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Abstract 

It is shown how a bare three-nucleon force is incorporated into the formahsm of the effective 
interaction approach for hyperspherical harmonics. As a practical example we calculate the ground 
state properties of and ^He using the Argonne V18 nucleon-nucleon potential and the Urbana 
IX three-nucleon force. A very good convergence of binding energies and matter radii is obtained. 
We also find a very good agreement of our results compared to other high precision calculations. 

PACS numbers: 21.45.+V, 21.30.Fe, 31.15.Ja, 21.10.Dr 
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I. INTRODUCTION 



It is well known that the binding energies of the three-body nuclei cannot be obtained 
with NN forces only and thus also three-nucleon (3N) forces have to be included. This makes 
realistic calculations for systems with more than two nucleons more difficult. Nonetheless 
such calculations have been carried out in three- and four-nucleon systems by various groups 
using different theoretical methods: Green Function Monte Carlo 111 (also for A > 41 as well 
as hyperspherical harmonics (HH) expansions and Faddeev techniques (see e.g. 
It would be very desirable to extend such an exact treatment of 3N forces to methods based 
on effective interaction approaches, since they are very suitable to treat also systems with 
A > 4. In case of the effective interaction approach within the No Core Shell Model (NCSM) 
such calculations have already been performed. Rather good results have been obtained ^, 
though without reaching a convergence as excellent as in other methods. For the effective 

r-i 

interaction hyperspherical harmonics (EIHII) technique |7| an inclusion of bare 3N forces 
had not yet been accomplished. In fact the aim of the present paper is to show how such a 3N 
force is incorporated in the EIHH formalism. A further aim is to check whether a comparably 
good convergence as with methods, which do not use effective interaction approaches, can be 
reached. To this end results of a calculation of three-nucleon groundstate properties using 
the AV18 NN potential '8] and the UrblX 3N force are presented. 

The paper is organized as follows. In section II the EIHII method is briefly outlined, 
while the incorporation of the 3N force in the EIHH method is described in section HI. 
Results for the three-nucleon groundstate properties are discussed in section IV. Details for 
the calculation of the 3N force matrix elements are given in the Appendix. 



II. BRIEF OUTLINE OF THE EIHH METHOD 

To apply an effective interaction method to an A-hodj Hamiltonian of the form 

i^I^l =Ho + V (1) 

one divides the Hilbert space of i/'"^^ into a model space P and a residual space Q with 

eigenprojectors P and Q of Hq. The Hamiltonian H^^'^ is then replaced by an effective 
model-space Hamiltonian 

ff[A]eff ^ pjj^p ^ py[A]effp ^2) 
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that by construction has the same energy levels as the low-lying states of H^^\ To find 
yi^e-ff ^ however, is as difficult as seeking the full-space solutions. In the EIHH method, as 
in the NCSM method, one approximates V^^^^^^ in such a way that it coincides with V for 
P — > 1, so that an enlargement of P leads to a convergence of the eigenenergies to the true 
values. 

Considering a Hamiltonian that contains besides the kinetic energy T only a two-body 
interaction vfj, 

H^^^=T + Y.-^, (3) 
one can rewrite H in the hyperspherical formalism as 

/7[^l = r, + T^(p) + V[^](p,fi^), (4) 

where 

V^^(,,a^^±^, T, = -J-A„ T,.(p) = i-M (5) 

are the bare two-body potential and the hyperradial/hypercentrifugal kinetic energies, re- 
spectively, m is the particle mass, Ka is the hyperspherical grand angular momentum op- 
erator and are the {3A — 4)-dimensional hyperangles, while Ap denotes the Laplace 
operator with respect to the hyperradial coordinate 



Y^VI with rf^ = ^ ^^^^ (^rl - f, j . (6) 

Here ffi are the N = A — 1 Jacobi vectors in the reversed order and fj are the position 
vectors of A particles. In the EIHII method one takes the so called adiabatic Hamiltonian 
as a starting point, 

n^^\p)^TM + v^^\p,nA), (7) 

and the unperturbed Hamiltonian Ho is chosen to be Tk{p), which has the HH basis functions 

aKj^] as eigenfunctions, where [Ka] stands for the appropriate set of quantum numbers (see 
). The P-space is defined as the complete set of HH basis functions with generalized 
angular momentum quantum number Ka < Kp and correspondingly the Q-space with 
Ka > Kp. 

For each hyperradius p an effective adiabatic Hamiltonian is constructed 

n^^^^"{p, Qa) = PTk{p)P + P V^^^^ff{p, Qa) p , (8) 



where \/[^]'^-^-^ is approximated by a sum of two-body terms 

A 

i<j 

Due to the use of antisymmetric wave functions one needs to calculate vff^^ only relative 
to one pair, since 

( v^*" > - { E > = > • (10) 

i<j 

For the determination of v'a^a-i) "-"^^ takes a p-dependent quasi two-body adiabatic 
Hamiltonian 

n^^\p;9N,m) = Tk{p) + v^^\^_,^{V2psm9N ■ m) , (11) 
where the unit vector fj^ and the hyperangle 6^ are defined by 

Vn = —, Vn = \ -ta), Vn = P sin On ■ (12) 

The Hamiltonian of Eg. (II ID is diagonalized on the A-body HH basis, which then allows to 
apply the Lee-Suzuki 10| similarity transformation to 'h}'^\p ; On, Vn)- This finally leads to 
an expression for a p-dependent two-body effective interaction v^'^^^^^{p ; 6^, Vn) (for details 
see j?!), and thus the A-body problem can be solved in the P-space with the following 
effective Hamiltonian 



i<j 

One repeats the procedure enlarging the P-s pac e up to a convergence of the low-lying 
energies of an A-body system. As shown in [ll| the EIHH method leads to a very fast 
convergence for the binding energies of few-nucleon systems with realistic NN potential 
models. 
In 



12| we could show that the convergence is even further improved if one goes beyond 
the adiabatic approximation taking into account also an effective hyperradial kinetic energy 
/\j'e//_ xhis leads to the following nonadiabatic two-body effective interaction 



For explicit expressions of IST^^^ we refer to 
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III. INCORPORATION OF THREE NUCLEON FORCE 



The inclusion of a bare three-body interaction v-j)^ leads to the following HH A-hodj 
Hamiltonian 

H^^=T, + TM + f:v^+j:vg. (15) 

i<j i<j<k 

While for the two-body interaction one needs to consider only the hyperspherical coordinates 
connected to the A-{A — 1) pair explicitly, for the three-body interaction one has to take into 
account additional coordinates. We were confronted with the same problem when deriving 
in 1^ a three-body effective interaction from a bare two-nucleon potential. Thus we can 
proceed here in the same manner. 

In order to focus on the interacting three-body subsystem we transform the (3^4 — 4) 
hyperangular coordinates Qa = {^2, ■ ■ ■ , ^71, ■ ■ • , Va-i) into a new set of hyperangles 

n3,A-3 = (0,, 9a-3, e3'^'^ 9l„ fj,,..., fjA-l) . (16) 

The new hyperangles reflect the splitting of the A-body system into 3- and {A — 3)-body 
subsystems. The hyperangles fls^A-s can be written as 

f^3,^-3 = (e[^'^-^(^EU[i-^i) (17) 

[31 

with the hyperangles of the interacting three-body subsystem {^int) ^-nd those of the residual 
system (f2i^7^') defined as 

^^t = (CI VA-2, VA-l) , ^It'^ = (02, . . . , ^A-3, , ^1, . . . , VA^S) ■ (18) 

The two new angles 9'^''^"'^^ and 6'|^'^, replacing 6a-2 and 6a-i, are given by the relations 
pPJ, ^ + v\-2 = P sin , pA-s = p cos Q^''^-'^ (19) 

and 

VA-i = pflt sin , Va-2 = & cos CI • (20) 

The new coordinates (pfnt^^fnt) forai a complete set for the three-body problem and thus 
matrix elements of 3N forces can be calculated for such a system. 

The calculation proceeds as follows. In a first step one solves the A-body Schrodinger 
equation with the Hamiltonian of Eq. and constructs the nonadiabatic two-body effec- 
tive interaction v^^^^^^^ . As next step one has two alternatives. One of them is to solve the 



Schrodinger equation with the Hamiltonian 

+ E4- (21) 

i<j i<j<k 

This is the effective Hamiltonian, which is used in the three-nucleon calculations presented 
in section IV. 

The other alternative (meaningful for A > 3 only) is to solve first the adiabatic quasi 
three-body Hamiltonian 

= TK{p) + V^'\pl„nl,). (22) 

The so defined V^^^ contains all three effective pair interactions fj^''^'^'^ and different from |l2l| 
also the bare 3N interaction v^^\. Then one can proceed further in complete analogy to |l2l |. 
i.e. along the same lines as for the construction of v^^"^^^^ applying the Lee-Suzuki similarity 
transformation in order to obtain v^^^^^^. Finally one solves the A-body Schrodinger equation 

H^^^^ff = T, + n^^^^ff c^T, + Tk+J2 ^Sk" (23) 

i<j<k 

taking into account that 

(V^[^1^^0^( E ^Sf ) = -^(-^ - - 2) ^ .,3j,,,^^^ A-l,A-2)). (24) 

i<j<k 

IV. DISCUSSION OF RESULTS 

In order to check the convergence behavior of the EIHH method in presence of 3N forces 
we choose as test cases the ground states of ^H and ^He. We use the AV18 NN potential ^ 
including the electromagnetic corrections of the NN force, but neglecting the isospin mixing. 
As 3N force we take the UrblX model jol. The calculation of the 3N force matrix elements 
is described in the Appendix. 

In Fig. 1 we show the binding energy of the triton as function of the grand angular 
momentum quantum number K. One observes a very good convergence pattern. In fact 
already with K = 8 one has only a difference of about 50 keV with the value of K = 20. 
This difference is further reduced to about 20 keV with K = 10, while with i^' = 12 one has 
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essentially the final result. In Fig. 2 we show the corresponding results for ^He. One finds a 
very similar picture as in Fig. 1, i.e. again an excellent convergence. In Fig. 3 we illustrate 
the results for the matter radii of and '^He. Also in these cases a very good convergence 
is evident with almost identical convergence patterns for '^H and ^He. 

In Table I we list our results for the triton binding energy and the probabilities for the 
different angular momentum components together with those of Table I in Q], where the 
inelastic longitudinal electron scattering form factors of the three-nucleon systems were 
calculated (for this comparison the electromagnetic part of AVI 8 is not included). There is 
a good agreement for the various probabilities and a small difference of about 40 keV for 
the binding energy. 

With an inclusion of the electromagnetic interaction of AV18 we obtain the results given in 
Table II, which are compared to those of The comparison of the various results reveals 
an excellent agreement for the various ground state properties obtained in the different 
calculations, e.g., the maximal difference for the binding energy is 6 keV for both three- 
body nuclei. Thus we may conclude that 3N forces can be included in the EIHH formalism 
with high precision. This is an important finding opening up the way to include 3N forces 
also for systems with A > 3 in this approach. 



Appendix: Multipole expansion and matrix elements of 3N forces 

We consider the Urbana and Tucson-Melbourne [3] versions of 3N forces. To calculate 
their matrix elements (ME) we first represent the 3N force operators as expansions over 
irreducible space tensors. We shall list these expansions in two different forms. 

The Urbana or Tucson- Melbourne 3N forces i;[23 given by a sum Ylcyc ^ijk taken over 
cyclic permutations of particles 1,2,3. Each of the three terms in the sum gives the same 
contribution to ME between states antisymmetric with respect to nucleon permutations. 
Each of the terms Vijk is also symmetric with respect to interchange of two of the three 
nucleons (in our notation they are nucleons i and j). When Jacobi vectors are constructed 
in the natural order the spatial components of basis states are chosen to be symmetric or 
antisymmetric with respect to the interchange of nucleons 1 and 2. Then it is convenient 
to use the V123 component of the force for the calculation of ME. The expression for this 
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component is of the form 



Vl23 = (^27r{Xl3, X23} + Bi23{zi3, Z23}) {n ■ T3, T2 ■ T3} 

+C27r[a;i3, X23] [n ■ T3, T2 ■ n] + Uo (Ti3T23f . (A.l) 

Here [...,...] and {. . . , . . .} mean commutator and anticommutator, respectively, while 

where Yij = Y{rij) and = T{rij) are the regularized Yukawa and tensor functions, 
respectively. As customary Sij is the tensor operator, and A2tt, C2-k, and Uq are the strength 
coustaMs (for the Urbaua models C.. = (1/4)^..). The numerieal values of A.^ a„d f/„ 
that correspond to the Urbana IX model are listed in j9[. The quantity Bijk = B{rik,rjk) is 
different from zero only in the Tucson- Melbourne model . 

In what follows we use the notation cTj and Xj for nucleon spin and isospin operators. We 
use the following notation for tensor products of irreducible tensor operators 

{Aa ® Bb)xi, = Yl {aOibp\Xfl)AaaBbp, {Aa ■Ba) = J2 AaaBa,-a{-lT, 

a+/3=fj, a 

(Aa^Bh)'^^ = J2 {aab/3\Xfi)[Aaa, Bhpl ® 5^)^^ = J] (aa6/?|A/i){A„„, S;,^}, 

and also the notation for irreducible spin tensors 

Use of spherical harmonics and their tensor products as irreducible space tensors gives 
the expansion of V123 in the form 

3 

^123 = Wo+ (F2(ri3) ■ WD + {Y2{r23) ■ WD + Yl (^A '(^3, r23) ■ WD , (A.3) 

A=0 

where Y^^{f is, f 23) = (^2(^13) ® ^2(^23))a^5 % stands for the unit relative vector of nucleons 
i and j {fij = ri — rj). The quantities Wq, W^^^, W^2^"^ W^^ are irreducible tensors with 
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respect to spin variables of the ranks 0, 2, and A, respectively. One has 



Wo 



f/o(Ti3T23)^ + (A2^yi3^23 + 5l23){n " ^3, f2 ■ f'i}{ai ■ (T3, ^2 ' 0^3} 
+C27r>l3^23[Tl ■ ^3, ^"2 " ^3] [(Ji ■ (T3, (72 " 0^3], 



247r 



[(A2^ri3r23 + 5i23){n ■ ^"3, T2 ■ t^}{t.^?,, 02 ■ 0^3} 



+ C'2,rTi3F23[n • r3, r2 • r3][Si|^, (72 ■ ^3] 

^(-1)^ [(A2.T13T23 + fii23){ri ■ r3, r2 ■ r3}(S23 ® S^g)^^ 
+C^2.Ti3T23[fi ■ r3, r2 • r3](S?3 ® S^3)^J , 



(A.4) 



and = (12)14^^^. The formulae above are obtained using the following expression for 

the tensor operator 



and the recoupling 



247r 



(F2(ri3) ■ S?3)(r2(r23) • ^Is) = E(-l)'(^A'(^i3,r23) • (S?3 ® S^3),) 

A=0 

(for A = 4 the quantity (S^3 ® T,l.^)x^ vanishes). 

The ME of the scalar products entering ()A.3jl are expressed in the standard way in terms 
of spatial reduced ME and spin-isospin ME reduced with respect to spin. The latter ME 
may be calculated directly writing down all the quantities depending on spin and isospin in 
terms of crj;o,± and rj.o,±. One may also employ the expansions over the tensors of Eq. ()A.2|) 
form with subsequent use of the ME flA.10|) . ()A.11|) . To obtain these expansions we perform 
recouplings, make use of the values of {o"3/i, o_^} and [0-3^, asi,] 3] and also perform simple 



summations of Clebsch-Gordan coefficients As a result we get the following expression 
for the 3N force 



V123 = 4y3(ri®r2)o{(A2.ri3V'23 + 5i23)v^S°2 



247r 



(1 + (12)) (A2^ri3F23 + Sl23)(>^2(ri3) " S?^) 



247r 



+5— (A2.T13T23 + 5i23) J2 \/3(2A+l) < 



A=0 



1 1 2 
1 1 2 
A A 



+ (1 + (12)) Ti3r23 {Y2ih,) ■ (sj^^g + Vssgg; 
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A=0 



1 1 2 
1 1 2 
A 1 A 



^ ^12,3 



+Uo{TuT, 
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(A.5) 



Another multipole expansion of the force is obtained when the quantities 



((rijOrij)A® {fik0rii)x') 



AM 



(A.6) 



are used as spatial tensors and, moreover, an expansion of spin-isospin operators over the 
tensors of Eq. ()A.2|) form is performed. In this case we assume use of Jacobi vectors 
constructed in the reversed order and for the calculation of ME we employ the component 
V231 of the force. We rewrite it in the following form: 

^231 = ^27r(2r2 ■ T3) {Ftt{^2 ' h2){^?. ' ri3)(ri2 ■ ri3) + Fyy{^2 ■ ^3) 
+Fty{^2 ■ ^12) (0^3 ■ ^12) + FYT{a2 ■ ri3)(CT3 • fi^)) 
-C2^(2ri ■ f2 X T^){Ftt{^2 ■ ^12)(o'3 ■ '^13)(o'l " X fi3) + Fyy(ai ■ (T2 X 0^3) 

+Fty{^2 ■ ri2){a-i ■ di X fi2) + Fyt^^s ■ ^13) ■ ^2 x hs)) 
+ 523i(2r2 ■ r3)18(a2 • ri2)(a^3 ■ ri3)(ri2 ■ hs) + UoiT^T^sf (A.7) 

with Ftt = I8T12T13 and 

Fyy = 2(ri2 - Ti2)(Fi3 - ri3) , Fty = 6Ti2(ri3 - T13) , Fyt = 6(Fi2 - Ti2)Ti3 . (A.8) 

After a recoupling of the various operators one obtains 
^231 = -2y3A2.(r2 ® r3)o{s°3 ■ [-V3{Ftt + 18^)X^,X^, - ^(3Fyy + Fty + Fyt)] 



A2^ 

231^ ^0 1 
^23J 



A2^ 

B 



+ E23 ■ [-\/?>{Ftt + 18^^^) -^23-^^23 + {FtyX22 + FytX^^)]^ 



A2n 

-4y3C2.(ri®(r2®r3)i)o<|s^3°i- v^[<; ^ ^ ^ \ {Fty + Fyt) - Fyy - 

1 1 I V 3 



^23,23 J 



+ ^23!l ■ [-^TT Ar^3!23l + ^23;i ' [-^TtX 



'(1,0)1 



23,23 
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+ ^23,1 





fill) 












1 1 2 








1 2 1 



23,23 



(1.2)2n 
23 J 



+Eg, ■ [FttXH^']} + U,{T,,T,,f. (A.9) 
For the spin operators of Eq. ()A.2|) one finds the following reduced matrix elements 

((Si; (S2; S3)S'23)5'i23| 1(0^2 ® 0^3)a||(si; (s2; S3)5'23)S'i23) 



5'23 'S'123 \ 



'5'l23 '^'23 A 



X 



7(2523 + 1)(2^^3 + 1)(2A + 1) < 



1 1 c 

2 2 "^23 

lie/ 

2 2 "^23 
1 1 A 



(A.IO) 



((Si; (s2;S3)S'23)5'i23||(ai (g) (^2 ® a3)A)A||(si; (S2; S3)5'23)5'[ 



123/ 



6 V(2^i23 + 1)(2S;23 + 1)(2A + 1) < 



I 5*23 5123 

1 01 01 

2 '-'23 '-'123 



X 



7(2523 + 1)(25^3 + 1)(2A + 1)< 



1 A A 

1 1 



lie 
22 ' 



23 



lie/ 
2 2 "^23 

1 1 A 



> . 



(A.ll) 



The isospin operators of Eqs. ()A.9|) and ()A.5|) lead to corresponding expressions for the 
reduced matrix elements. 

For the calculation of the spatial matrix elements one has to perform a six-dimensional 

Q 

integration. For this integration we use the same technique as first used in [2| and explicitly 
described in |16i]. For the Jacobi vectors (see section III) 



VA-l = pfnt sin VA-l , VA-2 = pfnt COS VA-2 



(A.12) 



one chooses a special reference frame in which 17^-1 turns into a vector (a-i = z and fjA-2 
turns into a vector Ca-2 that lies in the (x-z) plane. Then one replaces the angles 6*]^'^, r/A-i, 
and r)^_2 by 

X = cos{2e^^t) , y = fiA-i ■ flA-2 , (A.13) 
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and the three Euler angles to. This leads to the following volume element 

J3] \5 J^[3] \^fT~~3,Arr.A^ . .rif, . . — ( M ^[3] ^ 



df}A-idf]A-2 = {p\^tfdp\^t-Vl^^dxdfiA-idfiA-2 = {Pintfdp\^t-Vl - x^dxdyduj . (A.14) 
With this parametrization one obtains for the three-body HH functions 

YKtMi^) = i^%''{x)[Yt,mAVA-l)®Y,^m,{flA-2)]LM 

M' 

= E D^m'm{^)Y'/lm'{x, y) . (A.15) 

M' 

The five-dimensional integration in dxdyduj can be reduced to a two-dimensional one. In 

fact after some algebra one finds for a coordinate space operator T^A/u(pl^t, ^flt) T^Siuk A 
and projection p the following reduced matrix element 



MM'n 

xn' j dxVl^dyY^^Mi^. y)VxMSlt. y)>^'f'M'(^, v) ■ (A.16) 
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TABLE I: '^H ground state properties with AV18+UrbIX (no electromagnetic interaction, no isospin 
mixing) for binding energy BE and probabilities of total orbital angular momentum components 
in % 





this work 


CHH \3] 


BE [MeV] 


8.508 


8.47 


S-wave 


89.504 


89.55 


S'-wave 


1.061 


1.05 


P-wave 


0.135 


0.13 


D-wave 


9.300 


9.27 



TABLE IL Same as Table I but with inclusion of electromagnetic interaction (no isospin mixing) 





this work 


Bochum [4] 


Pisa [4] 


BE [MeV] 


8.468 


8.470 


8.474 


S-wave 


89.516 


89.512 


89.509 


S'-wave 


1.059 


1.051 


1.055 


P-wave 


0.135 


0.135 


0.135 


D-wave 


9.291 


9.302 


9.301 


TABLE IIL Same as 


Table II but for ^He 




this work 


Bochum [4] 


Pisa [4] 


BE [MeV] 


7.736 


7.738 


7.742 


S-wave 


89.385 


89.391 


89.378 


S'-wave 


1.246 


1.229 


1.242 


P-wave 


0.132 


0.132 


0.131 


D-wave 


9.238 


9.248 


9.249 
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FIG. 1: Triton binding energy with AV18+UrbIX (no isospin mixing) as function of grand angular 



mome] 




FIG. 2: Same as Fig. 1 but for ^He. 



15 



1.8 



1.6 



1 1 1 1 1 1 1 

• • ^He 

t — ^ ♦ ^ • • • i 




- - 
-^^.^-...-ff— ■ * - -■ i 




" / 

- / 
/ 

~/ 
It 

1 1 1 1 1 1 1 


1 


4 6 8 10 12 14 16 18 20 



K 



FIG. 3: and ^He matter radii with AV18+UrbIX (no isospin mixing) as function of grand 
angular momentum quantum number K. 
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